Mirror displacement energies and neutron skins. 
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A gross estimate of the neutron slcin [0.80(5) (A*' — Z)/A fm] is extracted from experimental proton 
radii, represented by a four parameter fit, and observed mirror displacement energies (CDE). The 
calculation of the latter relies on an accurately derived Coulomb energy and smooth averages of the 
charge symmetry breaking potentials constrained to state of the art values. The only free parameter 
is the neutron skin itself. The Nolen Schiffer anomaly is reduced to small deviations (rms=127 keV) 
that exhibit a secular trend. It is argued that with state of the art shell model calculations the 
anomaly should disappear. Ifighly accurate fits to proton radii emerge as a fringe benefit. 

PACS numbers: 21.10.Sf, 21.10.Ft, 21.10.Gv, 21.60.Cs 



Recent experiments H have considerably added to 
our knowledge of neutron radii, the most elusive of the 
fundamental properties of nuclear ground states. The 
two sets of measures are consistent with one another, 
and a recognizable pattern emerges Fig. 4], leading 
to an estimate for the neutron skin (v = neutrons, tt = 
protons, t = N — Z) 



Ar 



/(^ =-0.04(3) - 



1.01(15)^fm. (1) 



A third — totally different — experiment adds weight to 
this estimate: it deals with the sodium isotopes, lighter 
and far more exotic than the species studied in ij, |. 
Nonetheless, their A^^,^ behaviour is very much the same, 
as seen in Fig. |l|. 
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FIG. 1: Ar^„ from: [1, Fig. 4] (circles), the two approxima- 
tions of Ref. [3, Table I] (squares), and Eq. (uh 



On the theory side, we have an elegant analysis of 
A,._^^ @ I a-nd many mean-field calculations. Some of them 
are close to estimate (0), others not so close ||], though 
in the Na isotopes several Skyrme forces seem to give 
equally good results As was pointed out in 0, the 
neutron skins could be easily varied in such forces, with- 
out perturbing other observables, but the criteria to fix 
them were not obvious. The problem has been tackled 
head on recently P, 0| by constraining agreement with 



realistic calculations of the equation of state for neutron 
matter but no systematic study has appeared so far. 

The idea that started this investigation is that — 
assuming isospin conservation — a complete knowledge of 
proton radii would determine A^^^ . As the available ex- 
perimental data turn out to be consistent with a wide 
range of possible values of A^^^ , we resort to the Coulomb 
displacement energies (CDE) which are very sensitive to 
this quantity. The originality of our approach rests on 
the derivation, for all the observables, of smooth — Bethe 
Weizsacker type — forms that depend only on N and Z. 
The differences between these averages and the experi- 
mental values — i.e., the shell effects — will be found to 
exhibit secular behaviour that leads to high quality phe- 
nomenological fits for the radii. For the CDE, they lead 
to an optimistic assessment of the present status of a 
famous problem, the Nolen Schiffer anomaly [p^ . 

To relate neutron and proton radii we start by not- 
ing that rl = [E.<,(1/2 - il)(l/2 - ii)^] /A^ has the 
same form as the Coulomb energy (interchange rf^ by 
r-^) . It follows that (r^) = a{A,T) + (3{A,T)T, + 
7(v4, T)T^. This result relies on the same arguments that 
lead to the isobaric multiplet mass equation ||ll|, p. 302]. 
As (rj) is a functional of the state occupancies, a, (3 and 
7 may vary rapidly. To obtain a smooth form we as- 
sume some continuous occupancies, which make it possi- 
ble to obtain a four-parameter expression based only on 
dimensional considerations. Since we are interested only 
in mirror nuclei with — 2t, we propose 
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The overall A'^/'^ dependence is a general asymptotic re- 
sult for self-bound systems. The exp(g/^) correction ac- 
counts for the larger radii for small A. The terms in v 
and C, must be at most of the same order as po for large 
t/A and the only remaining uncertainty is in the scaling 
a in t/A"^. 

Obviously, we can derive a similar expression for pi, 
with t —t. Therefore v > represents a uniform 



2 



contraction of the two fluids, while C > imphes a tt- 
contraction and ^'-dilation, which gives the neutron skin: 



= A(C) = {py - p^) = 



(3) 



If (7 = 1 the volumes occupied by the two fluids may 
differ by a quantity of order A, and we have a "volume" 
skin. 

If (7 = 4/3 we have a "surface" skin, since the difference 
in volumes is at most of order A^^^. 

The volume option would be the correct one for strong 
attraction between like particles. In nuclei, the vn force 
is by far the strongest, and we adopt a surface skin. (A 
volume skin also leads to excellent radii, but it will be 
ruled out by the CDE) 

To determine the coefficients in Eq. (^) we shall fit 
the experimental values of charge radii (full references 
are given in reduced to point radii through the 

standard prescription p'^^P = [{rl)l^^j.g^ — 0.64]^/^ (used 
in H, full form in §). 

A fit is made to 82 experimental values for nuclei with 
iV or Z = 6, 14, 28, 50, 82, 126 (the classical closures 
which we call EI, for extruder- intruder), for which shell 
effects can be assumed to be minimal. The results are 
altogether remarkable: for 0.4 < C 1-2 the root mean 
square deviations (rmsd) are below 10 mf. We find nearly 
common parameters po — 0.943(2), g = 1.04(3), for dif- 
ferent (u , C) pairs, given in Table I for 0.6 < C ^ l-O- 

To decide which is the favored {v,C.), we compare 
63 experimental displacements, CDE=BE(Z->, A^<)- 
BE(Z<,7V>)— where Z> = 7V> = max(Z, iV)||||— with 
calculations. (We have followed where, out of 78 
values, those involving proton-unbound states are dis- 
carded). In the next paragraphs we define one by one 
the ingredients and then examine their effect column by 
column in Table I. 

The Coulomb potential Vc will be written in an os- 
cillator representation. To conform to standard use we 
introduce R-rr = (5/3) ^/^p^r- One of the many advantages 
of the oscillator basis is that one can separate trivially an 

adimensional two-body (l/r) operator through the first 
equality below 
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MeV, 



(4) 



The second equality provides the connection between 
CDE and radii. It follows from a famous estimate ||l^, 
Eq. (2-157)], which we apply separately to neutrons and 
protons (MeV units): 

hio^r 35.59 tiWi, 35.59 , . 

(5) 



(2Z)i/3 l^rl) ' (2iV)i/3 



{rl) 



The (l/r) operator will be evaluated in detail later, and 
it will be shown that Eq.(^ leads to a "smooth average 
of the diagonal monopole part" , 



MA{Z{Z ~1) - Z) 



MeV, 
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TABLE I: Root mean square deviations from observed CDE 
for three sets of coefficients in Eq. (^) with Aq = 16. Vc from 
Eq. (|l^). All energies in keV. See text. 
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identical — within a small exchange term — to the classical 
expression for the charged sphere, (i?^ is the charge, not 
the point radius, as explained after Eq. (|l^).) 

The column labelled (V^^) in Table I gives the rms 
deviations when only this term is included. The preferred 
value is at C < 0.6. 

The Bethe Weizsacker form for the charge symmetry 
breaking (CSB) potentials is 
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The parameters and Bs cannot be chosen arbitrarily. 
They are not known experimentally, but enormous theo- 
retical work has been devoted to CSB potentials, exten- 
sively described in Ref. ||l^, from which we take the Vb 
contribution to be « 100 keV around A = 16, and ~ 300 
keV in nuclear matter. Accordingly, we set = 16 and 
define a standard FJ* at B^ = 300 keV and B^ = 200 
VSn.- 



keV. Calling = Vq^ + Vg , the corresponding column 
shows that the preferred value moves to C = 0.8, with a 
slight gain in rmsd. 

In the last three columns By and Bg are allowed to vary 
simultaneously to give an idea of the range of plausible 
Vb parameters. (Note the advantage of introducing Aq in 
Eq. (|^): By — Bs equals the contribution a.t A — 16.) The 
constancy of the rmsd indicates the absolute need of con- 
straints on Vb ■ Conversely, the rapid variation of Vb con- 
strains the skin to a narrow range of plausible values for 
By and Bs- We propose[ p6| A ^,,, = 0.80(5) t/A, (nearly) 
compatible with Eq. (p]) p7| and (perfectly) compatible 
with the estimate for ^DSpj-, gj^ 

The results depend on the validity of Eq. (||) which 
was derived assuming isospin purity for the wave func- 
tions. The assumption is likely to hold near stability, but 
approaching the proton drip line it becomes more ques- 
tionable. This is why we followed Ref. in rejecting 
data from 15 pairs with proton- unbound Z^. If they are 
reincorporated, the numbers are telling: the preferred C 
moves slightly up at 0.85 but the rmsd doubles at 254 
keV. Some of the rejected pairs do fairly well, but for 
eight of them the error is above 500 keV. (In the re- 
stricted sample no error is larger than 300 keV). These 
findings are consistent with good T near stability and 
deteriorating away from it. At stability, the T « nuclei 
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raise the problem of a negative skin, small for A ^ 60 
and quickly dissappearing as T increases (see for exam- 
ple Our calculation should detect it as a slight 
deterioration for T = 1/2 with respect to T = 3/2 or 
higher. If anything, what is seen is the opposite. This 
observation does not contradict the results of Ref . jl^ ; it 
simply suggests that the hitherto neglected shell effects 
are stronger than those induced by T-violation. 

Let us examine now the smooth forms we have intro- 
duced. The radii enter only through the oscillator fre- 
quencies in Eq. (^), which define the basis in which we 
should work. For magic nuclei the basis reduces to a 
single state, and the observed radii can be assumed to 
be given by pjr- When moving away from closed shells, 
configuration mixing will increase the size of the system: 
moderately for nearly spherical nuclei, and radically for 
deformed ones. The crucial point is that p^^ is not sup- 
posed to produce good radii, but a good basis. To check 
whether this is the case we propose the following test: 
Try to fit all the available data through a shell corrected 
radius p^^ = + 1), fixing pj^ to the values we have ob- 
tained. The test is passed if the fit is good and nothing 
is gained by re-fitting p^. 
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FIG. 2: Comparison between experimental errors and differ- 
ences between observed and calculated radii = pn 
= 0.8, undistinguishable from other choices.) 



The definition of V is suggested by previous work on 
shell corrections We have 



(8) 



S stands for spherical, Q for deformed. The correction 
I? is a functional of the orbital occupancies in the EI 
(extruder-intruder) valence spaces. They consist in the 
orbits of harmonic oscillator shell p, except the largest 
(extruder), plus the largest (intruder) orbit from shell 
p+l. For protons, say, the EI degeneracy is D^^ = {p^^ + 
l){p^ + 2) + 2 (e.g., = 3 between Z = 28 and Z = 
50). The degeneracy of the non intruder orbits is D^tt = 
PT^{pTr + !)■ Now call z the number of valence protons, 
and define 5^ = z{Dt, - z)/Dl, = z{DrTv - z)/Dl; 



and similarly for neutrons. By construction, D vanishes 
at the EI closures. 

Fits of p^^ to 634 experimental values keeping the 
parameters previously obtained for p^^ — variation leaves 
them unchanged — lead again to remarkable results: with 
A = 5.6(2), p — 23(1) we obtain rmsd < 11 mf in all 
cases. Fig. S gives an idea of the quality of the results: 
most of the calculated points fall within (or very close to) 
the experimental error bars. The exceptions are the light 
nuclei (where halo orbits are important), and the region 
around the light Pt isotopes known for shape coexistence 
in the ground states. The test is passed: p^^ makes sense. 
It should be noted that a fit of Eq. (^) to all nuclei makes 
little sense. 

The Bethe Weiszaker form for Vb is obvious. 

The Coulomb potential is a different matter. The qual- 
ity and credibility of our approach rests on an exact treat- 
ment of V^j^. 

Following Ref. |^ we separate the Hamiltonian H = 
H„i + Hm- The monopole part Hm contains all terms 
in scalar products of fermion operators (a^ ■ Og) for sub- 
shells r and s. We consider the diagonal part involving 
the number operators = 0^-0^. The non-diagonal 
one — responsible for isospin impurities — will be ignored, 
as will the multipole Hm which contributes to shell ef- 
fects ]l||^|l|,||. _ 

Calling Viki'k' the matrix elements of l/r, the diagonal 
monopole part is 
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[J] — 2J + 1, Zk — number of protons in orbit k, where 
k = plj stands for the quantum numbers specifying a 
given harmonic oscillator (ho) orbit {p is the principal 
quantum number). The sum over the first k major shells 
containing r orbits can be reduced to a sum of fac- 
torable terms by diagonalizing the matrix ^{V^ifc}. This 
technique is extensively described in [^. As the high- 
est eigenvalue Er overwhelms all others, when combined 
with the corresponding eigenvector Uk and the single par- 
ticle counter term in Eq. (|^) it determines the largest 
contribution by far (zk is the number operator for orbit 
k) 




Zk Vkk 



(10) 



which amounts to a basically exact representation in 
Fock space: Er is 30 times bigger than the second 
biggest eigenvalue, and 100 times bigger than the third. 
Furthermore — as we shall see-Uk leads to an extremely 
coherent operator, while the other eigenvectors do not. 
We have treated the k = 8, t — 36 case, but the results 
are independent of the number of orbits. 

Fig. I shows the forms of Uk and Vkk = Vk- To deal 
with quantities of order 1, we have rescaled Uk =^ 



4 



\f{T)Uk and accordingly =^ t Er = .383. The 
(2j + l)-weighted averages over j-orbits within a major 
shell p are also shown {Up and V^). By referring the 
l{l + 1) term to its centroid we have to good approxima- 
tion 

Uk^Up- 0.01 [l{l + 1) - p{p + 3)/2]/(p + 3/2) (11) 

For Vfe the j dependence is more complicated. To leading 
order in a (p + x)~^ expansion the averages are Up ~ 
1.522(p + 1.4)-i/4, and Vp « 0.93(p + 3/2)-i/2. 

To extract a smooth form, we fit separately the two 
terms of Eq. ( [lO| ) at the p closures (i.e. replacing the sum 
over A; by a sum over p, with Zk = {p + l)(p + 2) = Dp, 
Uk = Up, and Vk = Vp) to forms a{Z — by. Exact values 
for U„, and Vp are kept. Then we insert the factor in 
Eq. m, to obtain (MeV units) 
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FIG. 3: The terms in Eq. (0). See text 



.858^1/6 [{Z - 1/2)2-1/6 _ ^1-1/6] 



(12) 



The numerical uncertainties of the fit 0.1%) allow the 
choice of the round numbers shown. This is the form used 
in the calculations. For clarity, in Eq. ^ it has been very 
slightly (and innocuously) simplified and the point radius 
parameter Rq = a/5/3po = 1.219 replaced by its charged 
value Rq = 1.226, so as to change the overall coefficient, 
which becomes the one for the classical charged sphere: 
Se^ZVSi?'^ « 0.864ZVi?'= « Q.IZ'^/A^'^. A pleasing 
result. 

The form of averaging we have used eliminates all shell 
effects, which show cleanly in the first panel of Fig. ^ as 
the difference between experimental and calculated val- 
ues. For clarity only t = 1 and 2 cases are kept. They are 
sufficient to indicate the secular nature of the deviations. 

Smooth filling can be replaced by ordely filling, which, 
using exact values for Uk and Vk in Eq. (^, should be 
close to a Hartree Fock (HE) result, since Coulomb ma- 
trix elements are quite unsensitive to details of the single 
particle wavefunctions, except in the case of halo orbits. 



The difference between smooth and HE approxima- 
tions is mostly due to the l{l + l) term in Eq. (11), whose 
influence can be understood by noting that above a closed 
shell with Z = Z^s it produces a single particle field 
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obtained from Eq. (|) and Eq. (0), using E^ = 0.383, 
and the coefficient in Up « 1.522(p-f 1.4)~i/'*, mentioned 



after Eq. (hd 
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FIG. 4; First panel: Experimental— calculated CDE for t = 1 
and 2. witht smooth occupancies [C, — 0.8). Second panel: 
As above for Hartree Fock occupancies. 



The HE displacements in the second panel of Fig. ^ 
indicate that the l{l + 1) effect leads to a dramatic im- 
provement for the A = 15-17 and A = 39-41 pairs, the 
two test cases most extensively studied in the literature 
([|0[ |l|, H and references therein). This gain is upset 
by a loss of regularity in the shell-effect patterns and the 
(apparent) need for a too strong Vb- If we keep V^ (and 
C = 0.8) the rmsd becomes large (330 keV). The figure 
corresponds to a fit (rmsd=150 keV) that lead to a huge 
change in Vb, at B., - = 425 - 180 = 245 keV. 

The situation is in all respects similar to that in 
Ref. 0, where a full HE description with the SKXcsb 
force also demands a too large Vb contribution (355 keV 
at ^ = 17). The agreement between ||lj, bottom Fig. 2c] 
and the second panel of Fig. ^ is quite good, to within an 
overall shift of some 100 keV for the latter A con- 
vincing indication that our approximation to a full HE 
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result is sound. However, it casts doubts on the very use 
of the HF approach: 

Near closed shells it is likely to be valid, and the improve- 
ment it brings about is very beneficial. Elsewhere, HF 
produces misleading patterns through shell effects that 
should be treated in conjunction with those of multipole 
and CSB origin. Unless this is done, it is clearly better 
to keep the smooth approximation. 

To conclude, we examine what can be said about the 
Nolen Schiffer anomaly (NSA), which we define as a sys- 
tematic failure to obtain good CDE without ad-hoc ad- 
justable parameters. The historical origin of the problem 
is clear from column Vq^ in Table I, showing that the 
Coulomb potential leads to large discrepancies for ac- 
ceptable values of the neutron skin It is also clear 
that the Vb potentials though small (of the order of the 
Coulomb exchange), are crucial in moving the skin to a 
reasonable range. As a consequence, the remaining dis- 
crepancies are the shell effects in the first panel of Fig. ^. 

As we have seen, for the classical test cases — at A = 
15, 17 and A = 39, 41 — the l{l + 1) term reduces the 
errors (to some 100 keV for the former and 50 keV for 
the latter). To achieve greater precision, good quality 
shell model calculations are necessary, a difficult task. 



not undertaken so far, at least for the CDE. Therefore, 
at present, we cannot decide whether there subsists an 
anomaly or not. 

For nuclei with several particles in a major shell, high 
quality configuration mixing is possible, and the task has 
been undertaken, not for the CDE, but for the differences 
of excitation energies between mirror yrast bands (CED, 
better called MED) Shell effects play a major role, 
and the CSB contribution is at least as important as the 
Coulomb one. The results achieve an accuracy of « 10 
keV on differences of up to 100 keV. This is the order of 
magnitude of the discrepancies we want to correct. Fur- 
thermore, the CED are unlikely to demand better control 
of the wavefunctions and the interactions than the MED 
do. As a consequence — though we cannot decide whether 
there remains a problem until the CED calculations are 
done — a fairly safe bet is that the Nolen Schiffer anomaly 
will disappear. 
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We shall use throughout Z> (^<) for "the proton (neu- 
tron) excess mirror partner". 

The error bars are probably too large. R. Machleidt con- 
siders By — 300 keV to be an upper limit (private com- 
munication) . 

It is often argued that measures are model dependent 
and sensitive to an "interaction" radius rather than to 
A,.„^ proper. Our answer to this objection is threefold. 
A) we are unaware of estimates indicating that uncer- 
tainties are larger than the — large — quoted errors. B) If 
anything, the effect would lead to smaller Ar^„ , as our 
calculations do indeed. C) We make no use whatsoever of 
the experimental results. The reader may decide for him- 
self whether he wants to compare them with our estimate 
or not. 

In [ ^ , bottom Fig. 2c] crosses and squares should be 
interchanged. Then it only takes good eyes to check the 
agreement. 

This statement of the problem is compact and accurate 
but non-standard. In the late 1970's some experiments 
appeared to be consistent with vanishing neutron skin, 
and it was pointed out that under this assumption the 
NSA would disappear |2^ . After this proposal turned out 
to be unrealistic, with few exceptions (such as [Q) little 
emphasis was put on the strong correlation between CDE 
and neutron skin. 



